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The new Hipparcos reduction 

van Leeuwen (2007) has produced a new reduction of the Hipparcos data that inchides a careful 
modeling of the satellite dynamics and eliminates the data correlations in the original catalog caused by attitude- 
modeling errors. As a result, the distances to the Pleiades and other clusters are now consistent with the values 
obtained with other methods. The new analysis also reduces the parallax uncertainties for the brightest stars 
by up to a factor of four. We use the new results to study the spatial distribution of massive stars in the solar 
neighborhood and to improve the distances to individual objects. 

The technique 

For a star with an observed parallax tTq and Gaussian uncertainty (7,^, its distance probability distri- 
bution is given by: 

1 - rTTn 



p{r\7ro) = Ar exp |^--^=— — j p{r), (1) 

where A is a normalization constant and p{r) is the underlying spatial distribution (which can be thought of as 
a Bayesian prior). Note that, in general: 

< p{rK) >^ 1/tt, (2) 

because of the conversion from parallax to its inverse, the availability of larger volumes of space behind I/ttq 
than in front of it, and the possibility of a non-constant p{r). This effect is known as the Lutz-Kclkcr bias (Lutz 
& Kelker 1973), who described it for the special case of constant p{r) (Figure 1). If p{r) is well-behaved (i.e. it 
has a cutoff for large r) and known, one can apply a Lutz-Kelker correction c such that: 



d = I/tTo + c, (3) 

where d is the median distance of p(r|7ro). Note that for a constant p(r), c is strictly infinite but can be made 
finite by introducing an artificial distance cutoff if (Ttt/tTo < 0.175. c is positive in most cases, with the exceptions 
arising when p{r) has a strong negative slope. 

In many astrophysical cases of interest p(r) is unknown a priori (e.g. one is measuring parallaxes 
precisely to obtain it). Under those circumstances, it may be possible to use the method described by Mafz 
Apellaniz (2001,2005). One starts by selecting a large homogeneous sample of objects with measured parallaxes 
and assumes a parameterized p{r), with the initial values for the parameters chosen from the literature. p(r|7ro) 
is then calculated for each star using Eqn.[l]and the results are used to build a new p(r) with a minimization 
algorithm that applies the binning recipe of Mafz Apellaniz & Ubeda (2005). The procedure is then iterated 
until the parameters converge. At the end of the process, one can check whether the choice for the p{r) 
parameterization was good and redo the full procedure with a new one if necessary. 
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Figure 1. p{r\TTo) for a constant p{r) and different values of ct^/tTo normalized to the value at r = I/tTq. The 
left panel has a linear horizontal scale and the right one a logarithmic one. Note how in the classical Lutz-Kelker 
case p(r|7ro) always diverges for large r. 
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Figure 2. Galactic vertical structure model and double semi-infinite cylinder used. marks the Sun position 
and the horizontal dashed line the Galactic plane. A sample p(r|7ro) (filled curve) and its projection on the 
vertical axis inside the considered volume (blue continuous line) are shown. Two additional sample projections 
ofp(r|7ro) are also shown (blue dashed lines). 
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Sample and parameterization 



At low Galactic latitudes b, the Hipparcos data are complete only for my < 7.9. Therefore, we can 
only study the massive stars in the solar vicinity. Using the data in the Hipparcos catalog itself, we selected 
the objects with an spectral classification of WR, O, or B0-B5, as well as all stars of luminosity type I. We 
also searched additional catalogs (e.g. van der Hucht 2001, Maiz Apellaniz et al. 2004) to check for possible 
omissions and analyzed the cases where the derived absolute magnitude was anomalously large in order to 
discard stars with erroneous spectral classifications (a total of 93). The final sample has 69 WR stars, 293 O 
stars without a WR companion, 871 non-0 stars with Imninosity class I, and 2770 non-supergiant B stars (a 
total of 4003 objects). 

Following Mafz Apellaniz (2001), we select a p{r) that depends only on the vertical Galactic coordinate 
z and that is composed of a self-gravitating, isothermal disk plus a Gaussian halo: 

p{z) = sech^ -f exp (i±gi) (4) 

or, equivalently: 

P{^) = - f ^ sech^ ( ^:±^] + ^ exp f '-±^) , (5) 

where pd,o and phfl ^ire the disk and halo volume number densities at z = 0, respectively; hd and are the disk 
scale height and halo half width, respectively; zq is the Sun's distance above the Galactic plane; a is the total 
surface number density; and / is the fraction of stars in the halo population. In order to minimize incompleteness 
effects due to extinction and distance when calculating the Galactic vertical structure parameters, we consider 
only the objects located in a double semi-infinite cylinder with a radius of 300 pc and a zone of avoidance given 
by l&l = 5° (Figure 2). 
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Figure 4. Massivc-star surface densities as a function of distance from the Sun vertical line. The black and red 
lines show the observed densities for all Galactic latitudes and for the region with |z| > 26.15 pc, respectively. 
The blue hne is the red line multiplied by a factor that accounts for the stars in the gap defined by \z\ < 26.15 pc. 

Results 

• The values of zq = 20.0 ± 2.9 pc and hd = 31.8 ± 1.6 pc (see Figure 3) are similar to the ones obtained 
by Maiz Apellaniz (2001) and other authors. 

• The results for the other three parameters are a — (1.91 ± 0.11) • 10^"^ stars pc^^, ~ 490 ± 170 pc, and 
/ = 0.039 ± 0.015. The low value of / indicates that the number of runaway stars in the sample cannot 
be too large. 

• There is little variation in the fitted results when selecting subsamplcs by spectral type or luminosity class. 

• There are large variations in the results when fitting subsamples along different Galactic quadrants. 
Specifically, Zq is larger for the third quadrant and smaller for the fourth one. Those results are caused by 
the presence of local disturbances in the distribution of massive stars (the Orion OBI and the Scorpius- 
Centaurus OB associations, see e.g. Eh'as et al. 2006). 

• There are also large variations in the massive-star surface density as a function of distance from the Sun 
along the Galactic plane (Figure 4). As previously noted by Mafz Apellaniz (2001), the Sun is placed at 
a local a minimum. A clear maximum exists at 100 — 160 pc, mostly due to the Scorpius-Centaurus OB 
association. At a distance of 200 pc, extinction starts to affect the observed population of massive stars 
near the Galactic equator and at 300 — 350 pc we reach the Hipparcos completeness limit even for higher 
latitudes. 

• The distance uncertainties for the nearest massive stars have been substantially reduced and the number 
of massive stars with significative Hipparcos distances has been increased. See Tables 1 and 2 and Figure 5 
for examples. 
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• The new Hipparcos distance to 7^ Vel is in excellent agreement with the recent values of SSG^^ pc derived 
with SUSI (North et al. 2007) and of 368^11 pc derived with VLTI (Millour et al 2007). Note also that 
the old Hipparcos value was less than 2 sigmas away from any of the current ones. It is often erroneously 
quoted that the old Hipparcos distance to 7^ Vel was 258^3]^ pc but that value docs not include the 
Lutz-Kelkcr correction c. 

• The new Hipparcos distance to 9^ Ori A is also in good agreement (1 sigma) with the VLB A distance to 
the Orion nebula of 414 ± 7 pc recently derived by Menten et al. (2007). 
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Figure 5. p{t\tto) for HD 188209 for three cases: [a] New values of ttq, cr^, and p{r), [h] old values of ttq, tT,r, 
and p{r), and [c] new values of tTq, and CTtt with constant p{r). In the three cases p(r|7ro) has been normalized 
to its peak value in the plotted range for clarity. 
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Old TT, 
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Old d 


Old c 




New TTo 


New d 


New c 






(mas) 




(pc) 


(pc) 




(mas) 


(pc) 


(pc) 


COph 


7 


12 ±0 


71 


144111 


4.0 


8 


92 ±0.21 


112+^ 


0.2 


6 Ori A 


3 


56 ±0 


83 


000 + 152 


42.3 


4 


72 ±0.58 


221+25 


9.5 


C Ori A 


3 


99 ±0 


79 




27.5 


4 


44 ±0.64 


2oq+43 


14.1 


15 Mon 


3 


19 ±0 


73 




200.3 


3 


54 ±0.50 




26.2 


CPup 


2 


33 ±0 


51 




115.3 


3 


00 ±0.10 


335+I1 


1.5 


72 Vel 


3 


88 ±0 


53 


278+^? 


20.5 


2 


99 ±0.32 


349+^^ 


14.4 


A Ori A 


3 


09 ±0 


78 


405+293 


80.9 


3 


03 ±0.55 




31.5 


fi Col 


2 


52 ±0 


55 




89.5 


2 


46 ±0.20 


412+32 


5.4 


^ Per 


1 


84 ±0 


70 




1369.8 


2 


61 ±0.52 


416tif 


33.3 


HD 149404 


1 


07 ±0 


89 


11801±g^t 


10866.3 


2 


40 ±0.36 


458+96 


40.9 


6*2 Ori A 


1 


72 ± 1 


00 


1867+g^f 


1285.3 


2 


11 ±0.42 




46.2 


10 Lac 


3 


08 ±0 


62 


353l^f 


28.1 


1 


88 ±0.22 




10.5 



Table 1. The twelve nearest O and WR stars with significant Hipparcos distances. 
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Name Old tTq Old d Old c New tTq New d New c 



(mas) (pc) (pc) (mas) (pc) (pc) 



Canopus 


10.43 


± 


0, 


.53 


97+^5 


0.9 


10, 


.56 


± 


0, 


.56 


96+g 


1.0 


Rigel 


4.22 


± 


0, 


.81 




17.2 


3, 


.78 


± 


0, 


.34 


267ti 


2.9 


Achernar 


22.68 


± 


0, 


.57 


utl 


0.1 


23 


.38 


± 


0, 


.57 


43^1 


0.1 


P Cen 


6.21 


± 


0, 


.56 




5.6 


8. 


.32 


± 


0, 


.50 


122t^ 


1.8 


a Cru 


10.17 


± 


0, 


.67 


motl 


1.8 


10, 


.07 


± 


0, 


.48 


100+^ 


0.9 


Betelgeuse 


7.63 


± 


1, 


.64 


176lif 


44.9 


6, 


.56 


± 


0, 


.83 




11.2 


Spica 


12.44 


± 


0, 


.86 


811^ 


1.0 


13 


.07 


± 


0, 


.69 


77tt 


0.6 


f3 Cru 


9.25 


± 


0, 


.61 


110+? 


1.9 


11, 


.70 


± 


0, 


.98 


88+« 


2.5 


Antares 


5.40 


± 


1, 


.68 


1027+?°°^ 


841.6 


5, 


.90 


± 


1, 


.00 




17.2 


Dencb 


1.01 


± 


0, 


.57 




12234.8 


2, 


.29 


± 


0, 


.32 


A75tll 


38.3 


e CMa A 


7.57 


± 


0, 


.57 




3.0 


8, 


.06 


± 


0, 


.14 




0.1 



Table 2. Results for massive stars with Vr < 1.5. 
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